Optical response in one dimensional Mott Insulators 
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We study the optical response of a Mott Hubbard system in the framework of the half-filled 
Extended Hubbard Model using the Density Matrix Renormalization Group (DMRG) method. We 
discuss the appearance of excitonic features inside the spectral gap as the system goes from the Spin 
Density Wave (SDW) to the Charge Density Wave (CDW) phase. 
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A detailed understanding of the optical response and 
charge gap in the one dimensional Mott insulators re- 
mains a challenge to existing theoretical methods. Re- 
newed interest in the subject stems from several experi- 
ments on materials such as SrCu02 and Ni halides [QJ^] 
with possible applications such as ultrafast switching in 
optoelectronic devices. The aim of this letter is to inves- 
tigate the effect of nearest neighbor Coulomb repulsion 
on the optical and Raman spectrum in these systems. 

A simplified model which has been used to describe 
the essential physics of these materials is the Extended 
Hubbard Model (EHM) defined as, 
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The first term corresponding to hopping between near- 
est neighbor sites and the second term to the onsite 
Coulomb repulsion provide the competition between itin- 
eracy and localization in the regular Hubbard model. 
The third term represents Coulomb repulsion between 
electrons occupying nearest neighbor sites. The Hamil- 
tonian as written above guarantees an insulating ground 
state with a filling of one electron per site. 

Although this model has been widely studied, many 
questions remain unanswered. The analytic approach at 
its best in Bethe-Ansatz provides an exact solution only 
for the V = case j|. The method has been used to 
obtain the energy spectrum and thermodynamics, but a 
reliable computation of dynamical quantities such as the 
optical response remains elusive. Non-perturbative ana- 
lytic studies of the dynamical response in these systems 
have largely been constrained to use the continuum limit 

8- 

Numerical methods such as exact diagonalization, al- 
though valuable in providing real frequency information, 
are limited to small system sizes [p|. Quantum Monte 
Carlo methods can treat large finite size clusters but an- 



alytic continuation from imaginary to real frequencies is 
an unreliable procedure. 

The renormalization group idea has helped deal with 
some of the toughest problems in physics characterized 
by large number of degrees of freedom playing an essen- 
tial role. The efficacy of the idea as a tool to compute 
experimentally relevant quantities lies in the ability to 
integrate out the non essential degrees of freedom. This 
has been brought to fruition with tremendous success 
in a numerical algorithm for low dimensional interact- 
ing systems known as the Density Matrix Renormaliza- 
tion Group (DMRG) ||. The numerical solution of fi- 
nite size systems on a computer is restricted by an expo- 
nentially increasing Hilbert space. The DMRG method 
works around this by a systematic truncation. It pre- 
scribes how to retain the most probable states of a sub- 
system required for an accurate description of a partic- 
ular set of states (usually, just the ground state) of the 
full system. The DMRG algorithm initially suited to deal 
only with ground state properties, was subsequently ex- 
tended to compute dynamical correlation functions . 
One of the approaches known as the "Lanczos vector 
method" constitutes choosing the particular set above as 
the ground state plus a set of Lanczos vectors. The vec- 
tors are chosen to approximate the reduced Hilbert space 
of excited states which connect to the ground state via 
the operator whose correlation function is desired. This 
method is very efficient in capturing low energy sharp fea- 
tures such as excitons in the optical spectrum; especially 
when bulk of the weight is in a single peak. Excitonic 
features in multi particle correlation functions have been 
observed in the EHM in previous studies §,0. 

The EHM at half-filling shows an interesting phase di- 
agram. In the weak coupling limit (U <C t) the system 
undergoes a second order transtion from a spin density 
wave (SDW) phase to a charge density wave (CDW) state 
as a function of increasing V, at V = U/2 + S(U). For 
intermediate values of U the SDW and CDW phases are 
separated by a narrow region with a bond charge density 
wave (BCDW) order. As one approaches strong coupling 
(U ^> t), the transition is again from an SDW to a CDW 
phase at V — U/2 + 5(U), but it is now first order. The 
small correction 6(U) is positive and approaches zero at 



1 



both the weak and strong coupling ends. The precise 
location of a tricritical point at the crossover between 
the first and second order transitions has been a subject 
of much investigation ]Tl|-[l5t and is complicated by the 
existence of the BCDW order. 

In our study we compute the optical response and lo- 
cal spectral function of the EHM in the strong coupling 
regime. We fix U at a realistic value of 12t. The first 
order transition between the SDW and CDW phases is 
manifest in the optical properties as well as in the ground 
state energy. The relevant values of V for SrCuC>2 and Ni 
halides are in the SDW phase. All the results presented 
in this letter were obtained from computations performed 
with finite size chains of N s — 50 sites with open bound- 
ary conditions using the Lanczos vector method. Study- 
ing other sizes (N s — 18, 34, 66) shows that results for 
N s = 50 are generic. We use the finite size version of the 
DMRG algorithm and choose m = 150 for the largest 
sizes. Selected runs performed with higher values of m 
did not introduce significant changes in the results. Typ- 
ical discarded weights were O(10~ 6 ). To validate our 
code we compare our results for the static and dynamic 
properties with exact diagonalization for short chains. 

We use the following definition for the response func- 
tion, 

XAB (u>) = 1 jf dte^+^iO | [At(t),B(0)] | 0} (2) 

The real part of the optical conductivity is defined 
through the imaginary part of the current-current re- 
sponse with e — > 0, 

o» = i;&M (3) 

where 

3 = ~it^2 ( c i+i<r c ^ ~ c ]* c j+i<?) ( 4 ) 

3,° 

is the paramagnetic current operator. The non-resonant 
Raman spectrum in a Mott- Hubbard system is given by 
the response function of the stress energy tensor Jl6| , 

In the case of an insulator \ TT is not the dominant 
contribution to the total Raman spectrum. But it is in- 
teresting for a comparison with Xjj because, j and r 
are respectively odd and even under parity conjugation, 
apart from an overall phase. 

In Fig.jj] we show the local spectral function for various 
values of V ranging from to 9t. In the SDW phase, the 
single particle gap (A s ) stays constant until a threshold 
value of V around 3t is reached and then starts reducing 
(cf. Fig.EI). Here and further in this letter when we refer 



to the gap in a correlation function we measure the po- 
sition of the lowest energy peak. We ignore the tail part 
which comes about due to the small finite broadening 
that is used to represent the Lanczos continued fraction. 
This implies that our values for the gaps are tight upper 
bounds to the actual ones. At the SDW-CDW transition 
the spectral gap reduces abruptly reaching a finite value, 
jumps up and then starts increasing again in the CDW 
phase. Note that in the CDW phase, the site on which 
we compute the spectral function is empty in the ground 
state. 
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FIG. 1. Local Spectral functions for different values of the 
coupling V. Note the abrupt change in particle-hole symme- 
try between the two phases. 

In Figj^ we report the current-current and non- 
resonant Raman response functions in the left and right 
columns respectively. A systematic change in the opti- 
cal response in the odd an d even (xtt) channels is 
discernable as we sweep through V across the SDW and 
into the CDW phase. For V = we see a broad feature 
centered at around U, in good agreement with a recent 
calculation for the standard Hubbard model within the 
DMRG approach |l7j] . Our method does not allow a res- 
olution of the tiny bump seen in the middle of the broad 
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optical absorption band. But, as V is increased we do no- 
tice the formation of a resonance which gradually gains in 
weight and shifts towards lower frequencies Jl(|[l8| . This 
constitutes a precursor of the excitonic feature that we 
describe further below. 
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FIG. 2. Optical and Raman response functions for different 
values of the coupling V. The vertical dotted line indicates 
the magnitude of the single-particle spectral gap in each case. 



For small values of V the optical (Ajj) and Raman 
(A TT ) gaps would be expected to coincide with the spec- 
tral gap, A s . We find them to be slightly larger because 
it is not possible to create fully non-interacting electron- 
hole pairs in a finite size chain. Beyond V ~ 1.5t, the 
optical gap falls below the spectral gap. The same hap- 
pens for the Raman gap around V ~ 3t. We define a 
quantity that we will call the excitonic weight (Wjj( TT ^) 
as the fraction of the weight in the optical (Raman) spec- 
trum below A s ; 



w _ it ' duJ Xj 3 {TT){u) 
vv n( TT ) r°° j / .\ 

J0 dUJ Xjj(TT)[U) 



(6) 



In Fig.|| we plot the excitonic weight as a function of 
V. In the case of both response functions it is seen that 
the excitonic weight is zero until the above mentioned 
crossing of gaps occurs (cf. Fig.||). As V is increased 
further, the excitonic weight starts appearing and a res- 
onance begins separating from the rest of the spectrum. 
When Wjj( TT ) reaches a maximum around V ~ 4i, the 
spectrum is dominated by a sharp excitonic feature car- 
rying most of the weight -86% (77%). This peak, well 
differentiated from the rest of the spectrum, is clearly 
located inside the single-particle spectral gap while the 
rest of the weight falls outside. The Lanczos method is 
rather well suited to describe these excitonic features, but 
it is not so good in capturing detail at the higher end of 
the spectrum. The optical bands inside the one-particle 
continuum tend to be shifted towards higher energies. 
At the same time the relative weight in the excitonic fea- 
tures is accurately represented, since we find that the sum 
rule for the optical conductivity in terms of kinetic en- 
ergy is obeyed with 1% accuracy or better (except very 
close to the transition). As V is increased further be- 
yond V ~ At the excitonic feature starts loosing weight 
and at the same time marches towards zero frequency. At 
the precise point of the SDW-CDW transition the exci- 
tonic mode reaches the lowest frequencies we can resolve 
(w ~ l/L). 
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FIG. 3. Excitonic weight as a function of the coupling V. 
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In a Mott insulator represented by the half-filled Hub- 
bard model, creation of an independent electron-hole 
pair (or a holon-antiholon pair in the Bethe Ansatz lan- 
guage) has a finite energy threshold; namely the spec- 
tral gap. This threshold is lowered in the presence of 
an attractive force by the binding energy of an electron- 
hole pair called an exciton. This attraction comes about 
due to the increased range of Coulomb repulsion in the 
EHM and is significantly absent in the standard Hub- 
bard model. As we increase V the energy gained in bind- 
ing electron-hole pairs keeps growing continuously and 
the excitonic feature moves closer and closer to zero fre- 
quency. At the same time the density of electron and hole 
states available for binding first increases, reaches a max- 
imum and then goes to zero at the SDW-CDW boundary. 
When the energy gained from binding the electron-hole 
pairs equals the energy cost of creating them across the 
single particle gap and the optical gap vanishes in both 
the odd (xjj) an d even (xtt) channels. 
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FIG. 4. Gaps for different values of the coupling V. 

Although our interest is in the strong coupling limit 
one can gain an understanding of the physics at hand 
by using the language of the continuum limit applied to 
the EHM p3|j ; an approach which is strictly valid only 
in the weak coupling case. Further using bosonization, 
the model can be split up into two Sine Gordon Models 
(SGM) for the charge and spin sectors respectively. The 
charge sector is the only one of interest as far as 
is concerned. On the other hand spin-charge separa- 
tion does not help a study of (Xtt) because the operator 
r involves both the charge and spin sectors. The ex- 
act solution of the SGM is available and the spectrum is 
governed by a coupling parameter usually called f3 [ ^9[ . 
For /3 2 < 47T the spectrum is built solely of kinks and 
antikinks. As (3 increases further, the SGM enters the 
attractive regime and kink-antikink bound states known 
as breathers are formed. These bound states in the SGM 
correspond to the excitons in the EHM that form as the 
value of V is increased. Given that we are interested in 
the strong coupling regime of the EHM, a quantitative 



comparison with the attractive regime of the SGM falls 
beyond the margin of applicability of the scaling limit. 

To conclude, we have shown that sharp excitonic fea- 
tures dominate the transport behavior in a particular 
regime of the Extended Hubbard Model at half-filling. 
These excitons are a direct consequence of the presence of 
non-local Coulomb interaction in this model. The Lanc- 
zos method combined with the DMRG approach is a pow- 
erful non-perturbative tool in computing the dynamical 
correlation functions of a non-trivial system. Several ma- 
terials such as SrCu0 2 , Ni halides and the (TMTTF) 2 X 
salts (where X is an inorganic monoanion) are believed 
to be good examples of Mott insulators |I|,p|,pO[ . Our 
numerical approach permits us to easily include other 
ingredients such as explicit dimerization and interchain 
hopping which are present in these materials in order to 
allow a quantitative comparison with experiments in the 
future. 
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